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Abstract. 

Recently, co mbining; radia l veloc ities from Keck/HIRES echelle spectra with published proper motion membership 
probabilities, ICote et al.l i2002h observed a sample of 21 stars, probable members of Palomar 13, a globular 
cluster in the Galactic halo. Their projected velocity dispersion ap — 2.2 it 0.4 kms~^ gives a mass-to-light ratio 
M/Lv = 40tl7, about one order of magnitude larger than the usual estimate for globular clusters. 
We present here radial velocities measured from three different CCD frames of commissioning observations ob- 
tained with the new ESO/VLT instrument FLAMES (Fibre Large Array Multi Element Spectrograph). From 
these data, now publicly available, we measure the homogeneous radial velocities of eight probable members of 
this globular cluster. A new projected velocity dispersion ap = 0.6-0.9 ± 0.3 kms~^ implies Palomar 13 mass-to- 
light ratio M/Lv = 3-7, similar to the usual value for globular clusters. We discuss briefly the two most obvious 
reasons for the previous unusual mass-to-light ratio finding: binaries, now clearly detected, and more homogeneous 
data from the multi-fibre FLAMES spectrograph. 
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^ 1. Introduction 



All the dynamical studies of nearby globular clusters have 
established that these dynamical systems contain no dark 
matter, apart from the expected stellar remnants such as 
white dwarves and neutron st ars (e.g., IPrvor fc MevlanI 
1X9931: iMevlan k Heggid llQQTj) . Consequently, globular 
clusters may be the most massive stellar systems in which 
no non-baryonic dark matter is dynamically detected, 
while dynamical evidence for non-baryonic dark matter 
seems to be present in most galaxies, from the faintest 
dwarf spheroidals (dSphs) to the brightest cDs galaxies, 
and clusters of galaxies as well. 

Some of the local dSphs, around the Galaxy and M31, 
have integrated absolute luminosities similar or fainter 
than those of the brightest Galactic globular clusters. 
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Since there is evidence that some Galactic dSph galax- 
ies are dark-matter-dominated, it is therefore reasonable 
to check if some globular clusters do present dynamical 
evidence for non-baryonic dark matter. Such a possibil- 
ity may be emphasized by the current predictions of Cold 
Dark Matter (CDM) numerical simulations of galaxy for- 
mation, in which the number of low-mass dark-matter sub- 
structures orbiting the halo of massive galaxies largely 
exceeds the number of dwarf galaxies ob served in the 
halos of both our Gala xy and M31 (e.g., iKlvpin et alJ 
I1999I: iMoore et all Il999. 2001). However, recently im- 
proved CDM models may c orrectly predi ct the observed 
number of satellite galaxies l)Binnevll2003j) . 

Could some of the globular clusters in the outer parts 
of the Galactic halo be such dark-matter substructure ? 
These remote stellar systems have so far been poorly stud- 
ied because of the difficulties in the acquisition of high- 
quality radial velocities and proper motions, direct con- 
sequences of the faintness and sparsity of these distant 
stellar systems. They are nevertheless important probes of 
the formation and evolution of the Galaxy, as their ages 
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and metallicities provide direct constraints on the dura- 
tion of halo formation process and on the time-scale for 
Galactic chemical enrichment, while the shape and extent 
of the Galactic dark halo are constrained by their orbital 
properties. 

In 1998, a program was started at Californian Institute 
of Technology to study the internal dynamics of seven 
distant halo globular clusters using the High Resolution 
Echelle Spectrometer (HIRES) at the W. M. Keck 
Observatory. The aim was the first direct measurements of 
the velocity dispersions and mass-to-light ratios for these 
clusters. Six clusters in this sample exhibited velocity dis- 
persions (Tp ~ 1 kms~^, translating into mass-to-light ra- 
tio values typical of globular clusters M / Ly ~ 3 (all mass- 
to-light ratios quoted in this paper are in solar units). 
Only one cluster, the halo globular cluster Palomar 13, 
dis played a velocity d ispersion larger than expected. 

bote et al.l ll2002l also referenced below as the Keck 
study) presented a careful analysis, combining radial ve- 
locities from Keck/HIRES echelle spectra wi th published 
prope r motion membership probabilities from lSiegel et all 
l|200l[) . They obtained a sample of 21 stars, probable mem- 
bers of Palomar 13. Their projected, intrinsic velocity dis- 
persion of (Tp = 2.2 ± 0.4 kms~^ implied a mass-to-light 
ratio = 401^7 : about one order of magnitude larger 

than the usual value for globular clusters. Cote et al. dis- 
cussed at length all possible reasons for such an unusual 
result: (i) some velocity "jitter" among the red giants, (ii) 
a few binary stars, (iii) a non-standard mass function, (iv) 
process of dissolving into the Galactic halo through catas- 
trophic tidal heating during a recent perigalacticon pas- 
sage, or (v) the presence of a massive non-baryonic dark 
matter halo. 

It is worth emphasizing that, in Cote et al., careful de- 
termination of the error bars made the usual mass-to-light 
ratio values for globular clusters Ai / Ly 3 at about two 
sigmas from the Ai/Ly value obtained for Palomar 13. 

Because of this marginally significant and puzzling re- 
sult from Cote et al., some more spectroscopic data of 
stars in the field of Palomar 13 were acquired during the 
commissioning of the ESO/VLT instrument FLAMES. We 
present hereafter new high-quality and homogeneous ra- 
dial velocities for 46 stars, 9 of them being members of 
Palomar 13, which provide new velocity dispersion and 
mass-to-light ratio values for this globular cluster. The re- 
maining of this paper is as follows: Sect. 2 presents the 
observation and the data reduction. Sect. 3 discusses the 
membership of the stars. Sect. 4 gives the new velocity dis- 
persion and corresponding mass-to-light ratio, and Sect. 5 
discusses the plausible reasons for the difference between 
the present results and those obtained by Cote et al. 

2. Observations and data reduction 

2.1. Observations 

All the new data presented in this paper were acquired be- 
tween August 28 and September 3 2002, during the com- 



missioning of the ESO/VLT instrument FLAMES (Fibre 
Large Array Multi Element Spectrograph), at Paranal, 
Chile. For more information on this instrument, see the 
ESO we b page httv://www.eso .orq and the recen t publica- 
tions bv lPasauini et all ll2002l) and lRover ei~an 0002). In 
order to clarify some possibly confusing uses of denomina- 
tions, it is worth mentioning that FLAMES is the name 
of the instrument, while GIRAFFE is the name of the 
spectrograph within FLAMES. 

We were involved neither in the selection of the stars 
nor in the preparation and acquisition of these observa- 
tions. Recently, these data were made publicly available 
from the ESO/VLT archives site^. 

The Palomar 13 observations are summarized in 
Table ^ Three exposures, each of one hour integration, 
were acquired in early September 2002. Calibration files 
were acquired a few days earlier. All data were taken in 
the high-resolution spectrograph setup HR-9 with spec- 
tral domain spanning from 514.3 to 535.4 nm and resolu- 
tion X/SX 20 000. Both positioning plates of MEDUSA 
were used with the 133 fibres assigned as follows: 55 fibres 
to program stars, 14 to sky measurements, 5 to simulta- 
neous wavelength calibrations, while 59 fibres remained 
unused. Average signal-to-noise ratio (SNR) for star with 

= 18 is given as indication of spectra exposure. Note 
that expected accuracy for HR-9 setup a nd SNR = 10 is 
cry, = 0.20 kms-i (see iRover et al.ll2002j) . 

2.2. Data reduction and calibration 

For each positioning plate, two calibration frames were 
used. The first one - the flat field - gave the master lo- 
calization and analytical model of the PSF perpendicu- 
lar to the spectral dispersion direction, while the second 
one - ThAr wavelength calibration frame - provided two 
pieces necessary for the master wavelength solution: the 
2-D optical solution and the 2-D Chebyshev polynomial 
correction. For the HR-9 setup, approximately 6 000 line 
positions were used to adjust the 30 parameters of the so- 
lution. Finally, the rebinned extracted ThAr spectra were 
cross-correlated with the ThAr mask and the measured 
mismatch between spectra was used to update the slit ge- 
ometry model. 

The spectra were extracted using the standard 
Python version of BLDRS - Baseline Data Reduction 
Software (girbldrs-1.09 and girbldrs-pipe-1.05) available 
from http: //qirbld rs. sourceforqe.net . B asic description of 
BLDRS is given in lBlecha etaP l)2000() . 

The most important calibration feature of GIRAFFE 
is the presence of five Simultaneous Wavelength 
Calibration (SIMCAL) ThAr spectra in each exposure. 
The SIMCAL spectra are optimally exposed and regu- 
larly spaced over the detector (spectra #1, 32, 63, 94, and 
125 out of 133 in the present setup). They are used to 
accurately adjust the master localization and the master 
wavelength calibration. 



^ ^http:/ /www. eso.org/science/flarnes_corn.rn/FLAMES_cornm_ PAL_lS.htrn\ 
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Table 1. Log-book of the observations from FLAMES commissioning data archives: Palomar 13 observations and 
corresponding calibrations collected in HR-9 setup (514.3-535.4 nm). SNR is the average signal-to-noise ratio for star 
of y- magnitude =18. 



# 


Archives file name 


Date of 


Time 


Instrument 


Exp. time 


SNR 


Remark 






observations 


[UT] 


configuration 


[s] 


at V = 18 




1 


FLAMES_GIRAF_OBS244_0015 


2002-09-01 


05:29 


Medusa2 


3600 


10 


good 


2 


FLAMES.GIRAF_OBS244_0016 


2002-09-01 


06:45 


Medusal 


3600 


10 


acceptable 


3 


FLAMES_GIRAF_OBS246_0003 


2002-09-03 


03:41 


Medusa2 


3600 


7 


poor 


4 


FLAMES_GIRAF_FLAT240_0008 


2002-08-28 


19:51 


Medusal 


143 


flat 




5 


FLAMES_GIRAF_WAVE240.0011 


2002-08-28 


20:10 


Medusal 


484 


ThAr 




6 


FLAMES.GIRAF_FLAT240.0020 


2002-08-28 


23:53 


Medusa2 


143 


fiat 




7 


FLAMES.GIRAF.WAVE238.0018 


2002-08-26 


19:33 


Medusa2 


482 


ThAr 





The raw images were processed through the follow- 
ing steps: bias subtraction, localization adjustment, opti- 
mal extraction, rebinning to linear wavelength space with 
step of 0.005 nm using the master wavelength solution, 
wavelength solution translation using the SIMCAL spec- 
tra cross-correlation with ThAr mask and the final rebin- 
ning using the translated wavelength solution. Note that 
we did not subtract dark (negligible effects) and we did not 
flat-field (irrelevant for the cross-correlation), the flat-field 
frame being only used to derive master localization and 
PSF model. It is worth mentioning that the final trans- 
lation of the master wavelength solution, which is below 
0.001 nm for all three science CCD frames, indicates an 
excellent instrument stability, at least over the period of 
one week spanning the above acquisition of data. 

2.3. Radial velocities 

Radial velocities 14 were then measured through the 
cross-correlation of the above spectra with the stan- 
dard GIRAFFE mask for the IIR-9 setup. It is a binary 
CORAVEL-type mask for the stellar spectral type FO, 
with variable line- width ranging from 0.005 to 0.06 nm and 
a total number of 98 spectral lines. In this paper, the new 
FLAMES data are labeled with F for FLAMES while the 
data from Cote et al. are labeled with K ior Keck. Table|21 
summarizes our results for all observed stars. Columns 
(1) and (2) give the star identification number and its 
cluster membership probability, first from this FLAMES 
study and second from the Keck st udy, wheneve r avail- 
able. All our stellar IDs come from ISiegel et~a 
Columns (3) through (6) give the equatorial coordinates 
a and d, the distance R of the star to the cluster cen- 
ter (aj2ooo = 23''06"M 4!480, S.imm = 12° 46' 19^.'20) and 
the V magnitude from ISieeel et al.1 l)200l[l . Columns (7) 
through (9) give the radial velocities, whenever available, 
from each of our three FLAMES epochs. Column (10) 
gives for each star the number n of valid measurements 
available from our three epochs, (11) indicates the cross- 
correlation-peak quality Q (1 for best and 9 for worst), 
where the first number represents the certainty of the 
peak identification and the second the quality of the fit, 
and (12) gives the standard deviation of the two or three 



measurements, whenever applicable. Column (13) gives 
the FLAMES mean radial velocity from our n measure- 
ments and the estimated error (cry^ or individual error) 
and Column (13) is the Keck mean radial velocity from 
Cote et al. for comparison. 

The raw measurements of 55 stars were cleaned by 
removing all measurements with unclear peak identifica- 
tion or very bad quality of the fit. From originally 55 ob- 
jects, 46 remain with at least one valid radial velocity 
measurement. We get 4 objects with all three measure- 
ments (^26), 30 objects with two measurements and 12 
objects with one measurement (6 from first plate, 5 from 
second and 1 from third plate). 

Ignoring the three Vr variable stars #36 #38 and #156 
(see discussion below), we have an immediate Vr quality 
check for the 27 -I- 4 = 31 stars with at least two measure- 
ments: the mean of their individual standard deviations is 
0.54 kms-i. 

We also include in Table [3 the 13 stars with one mea- 
surement only but having acceptable cross-correlation- 
peak quality Q (Column 11). Note that measurements 
with Q values higher than 2-3 should be considered with 
caution. 

Figures n and 13 give the finding charts of all 46 stars 
measured with FLAMES. In Table [21 all stars are sorted 
by increasing values of their mean radial velocities listed 
in Column (13). In this table, the intermediate horizontal 
line isolates the 11 stars (at the bottom of Table with 
very similar radial velocities: they may be considered as 
cluster members when using the radial-velocity member- 
ship criterion only (see below). 

It may be worth mentioning that all radial velocities 
listed in Columns (7), (8), and (9) are raw velocities, 
obtained directly from the cross-correlation process, cor- 
rected only for the solar system and earth velocities. No 
other attempt is made to correct for any effect due to var- 
ious dependence on stellar spectral types (the same cross- 
correlation mask was used for all spectra) or any effect due 
to the limited spectral range. The undergoing analysis of 
other FLAMES data sets not related to Palomar 13 in- 
dicates that the random error due to the mask mismatch 
could be a few 0.1 kms^^ while the systematic shifts due 
to the limited wavelength coverage do not exceed a similar 
value (a few 0.1 kms^^) while using the HR-9 setup. 
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Fig. 1 . Digital Sky Survey image with all 46 stars listed in Tabled Stars are numbered according to ISiegel et alJ 
i|2nnih . Member stars are enclosed by boxes and other program stars indicated with crosshairs. The entire FLAMES 
field (diameter 26') is shown. See Fig. [21 for a zoom into the cluster center. 



3. Memberships 

3.1. Membership from radial velocities 

A first glance at Column (13) of Table |21 reveals immedi- 
ately a clump of 11 stars whose radial velocities accumu- 
late around a velocity close to the systemic radial velocity 
of Palomar 13, namely 24.1 ± 0.5 kms~^, as measured by 
Cote et al. This clustering of radial velocities is significant: 
the Vr of all the other stars in Table El are away from the 
above systemic radial velocity by more than 10 times the 
standard deviation of the of these 11 stars. 

These 11 stars are potential cluster members, but be- 
fore computing any new estimates of the systemic ra- 
dial velocity and projected intrinsic velocity dispersion of 
Palomar 13, we first compare our measurements with the 
Keck radial velocities and check for possible non-members 
using proper motion and photometric data. 

3.2. Agreement with Keck radial velocities 

The last two columns in Table (21 provide, whenever avail- 
able, i.e. for 13 stars, the radial velocities from both 
FLAMES (K-FLAMES) and Keck (FrKeck) studies. The 
global agreement is excellent. From the 13 stars, 9 lie 



within the limits |KFLAMES- KKeck] < 1.2 kms"! with 
an average difference of 0.29 kms~^. 

From the six stars with both FLAMES and Keck mea- 
surements, which are clearly non-member (above the in- 
termediate horizontal line in Table only one, with 
FLAMES ID #156 (Keck ID #23), has a radial velocity 
difference between two FLAMES epochs which is larger 
than the above 1.2 kms~^ limit. This star, clearly de- 
tected as variable through our internal consistency check, 
has a mean KFLAMES = 7.79±2.11 kms"^ (mean of two 
measurements) and a V^Keck = 18.57 ± 0.82 kms~^ (one 
measurement), with \AVr\ = 12.27 kms^^. In both stud- 
ies, this star is not considered as a cluster member, nei- 
ther from both proper motion nor from Color-Magnitude- 
Diagram (CMD) criteria. We consider this star as variable 
according to radial velocity (either pulsating or binary). 

The situation is slightly different for the seven poten- 
tial cluster members (below the intermediate horizontal 
line in Table O having both FLAMES and Keck measure- 
ments. 

- The star with FLAMES ID #66 (Keck ID #41) has 
a V^FLAMES = 26.69 ± 0.5 kms^^ and a V^Keck = 
19.24±0.47kms"\ with |AK| = 7.45 kms'^ In both 
studies, this star is considered as a cluster member. It 
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Fig. 2. Detailed finding chart of the central field: y-band ima ge of Palomar 13 taken with the Low-Resolution Imaging 
Spectrometer on the Keck II telescope from lCote et alJ (2003)- Member stars are enclosed by boxes and other program 
stars indicated with crosses. 



is interesting to note that the above Keck value is the 
mean of three individual measurements, spread over 
more than a year, with KKeck = 18.91 ± 0.93 kms~^, 
18.37±1.84 kms-\ and 19.43±0.54 kms"\ in chrono- 
logical order. They do not show any significant varia- 
tion. 

There are two more stars with marginal disagreements: 
first, stars with FLAMES ID #37 (Keck ID #72) has 



a FLAMES = 26.32 ± 0.2 kms^^ and a FrKeck = 
28.79 ±0.27 kms"\ with |AVr| = 2.47 kms"^ sec- 
ond, star with FLAMES ID #88 (Keck ID #36) has 
a INFLAMES = 27.97 ± 0.5 kms^^ and a V^Keck = 
25.29 ±0.89 kms"\ with |AK | = 2.68 kms^^. 

Arc these three stars pulsating or binaries ? Two of 
them have only one FLAMES measurement, while the 
third star, with FLAMES ID #37 has two FLAMES V, 
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Table 2. Barycentric radial velocities for 46 stars in the field centered on Palomar 13. The extra horizontal line isolates 
a clump of 11 stars with radial velocities very close to the systemic radial velocity of the star cluster, potentially 
members of the cluster. 
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measurements differing by 0.54 kms^^. We discuss below 3.3. One RR-Lyrae Star in the cluster 
the influence of these stars on the velocity dispersion es- 
timate of the cluster. 

The star with FLAMES ID #38 has two measurements, 
with FLAMES = 35.87 ± 0.70 kms^^ and 6.92 ± 
1.10 kms"\ with |AVr| = 28 95 kms~^ This is a known 
RR-Lyrae variable (##3 in ICiatti et~a'lj.l96,5.) . however 
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we ignore it in the remaining of this study in spite of the 
fact that it is considered as a cluster member. 

3.4. Membership from stellar proper motions and 
photometry 

Table 13 gives for the clump of 11 stars (at the bottom of 
Table [SJ, all of them potential cluster mem bers through 
their V^, the following information (from ISiegel et alJ 
l200l[) . Columns (1), (2), and (3) give the ID number and 
the equatorial coordinates a and S. Columns (4), (5), and 
(6) provide the proper motion values fi{a), /i(<5), and their 
related membership probability in %. Column (7) gives the 
FLAMES radial velocity. Columns (8), (9) and (10) dis- 
play the U,B,V photometry, while Column (11) indicates 
the position of the star on the CMD: RRL for RR-Lyrae 
pulsating star, RGB for Red-Giant Branch, SGB for Sub- 
Giant Branch, BS for Blue Straggler, and off-r for being 
off to the right of the cluster sequences. 

The me r nbers hip probabilities as determined by 
ISiegel et al.l l)200l|) are based on the proper motion rel- 
ative to the average proper motion of the predetermined 
cluster members. The probability is computed using, for 
each individual candidate star, the ratio between the rela- 
tive proper motion /i and the measurement error a^. The 
process is iterated till a clean sample of cluster members 
covering the area of study is obtained. 

The proper motion membership probabilities in 
Table O give 8 stars out of 11 as clear members of the 
cluster, while the 2 stars with FLAMES ID #24 and #127 
have zero probability of being members. The RR-Lyrae, 
FLAMES ID #38, with a low membership probabihty of 
48 % is ignored anyway. 

The CMD position in Table|3|gives 8 stars out of 11 as 
clear members of the cluster, corresponding to the clus- 
ter RGB, SGB and BS sequences. Again, the same two 
stars with FLAMES ID #24 and #127 are off the cluster 
sequences, to the right of the CMD, although less signif- 
icantly for the former than the latter. This confirms the 
selection based on the proper motion. 

The star with FLAMES ID #79 (not measured in 
the Keck study) has a CMD position indicating that it 
may be a blue straggler; given its T^FLAMES = 25.76 ± 
0.65 kms~^ and its proper motion membership probabil- 
ity of 65 %, we consider it as a cluster member. 

3.5. Eight members of Palomar 13 7 

The above discussion illustrates the danger of considering 
only radial velocities in order to select cluster members. 
The two stars with FLAMES ID #24 and #127 (the lat- 
ter not measured in the Keck study) with VrFLAMES = 
26.44±0.42 kms^i and 25.39±0.21 kms'^ respectively, 
have Vr values very close to the cluster systemic radial 
velocity. However, both proper motions probabilities and 
CMD positions make us consider them as non-members of 
the cluster. 



From the examples above, it becomes obvious that to 
declare a star as a member or not of Palomar 13 represents 
a difficult task, which needs the compulsory simultaneous 
use of at least the three sorting processes available so far, 
viz radial velocities, proper motions, and positions in the 
CMD. 

We are left with a sample of eight stars probable mem- 
ber of Palomar 13. 

4. Velocity dispersion and mass-to-light ratio 

Given the limited size of our sample which docs not permit 
any attempt to discuss any concentric variation of any 
kinematic quantity, we focus here on the best estimate 
of the unbiased mean radial velocity dispersion for the 
whole sample. The crucial element is the subtraction of the 
random error of our sample and the error estimate of our 
velocity dispersion. We compute the velocity dispersion 
and the mass-to-light ratio in four different ways: (i) for 
the eight stars with at least one FLAMES Vr value, (ii) 
for the five stars with two FLAMES values (measured 
twice), and then for the six stars which are measured in 
both Keck and FLAMES studies, (iii) for the six FLAMES 
Vi values and (iv) for the six Keck Vr values. 

The average empirical variance of radial velocity of 
31 stars with at least two measurements is varQ = 
0.29 km^s~^. This is the most realistic estimate of the 
measurement error applicable to all measurements. The 
stars with two measurements are taken with var = 
O.bvarQ while stars with only one measurement with 
var — varQ. We compute the unbiased systemic cluster 
velocity (Vr), the intrinsic projected cluster velocity dis- 
persion (Tp and asso ciated errors according to the method 
described in lPrvor fc Mevlan (1993,) . which were also used 
in the Keck study (Cote et al.). 

4.1. Sample of eight stars with at least one FLAMES 
Vr value 

The systemic velocity of Palomar 13 from our eight mem- 
ber stars is (K-s) = 25.91 ± 0.36 kms~^ with a projected 
velocity dispersion Cps = 0.92 ± 0.29 kms^^. This (jpg 
value should be considered as an upper limit to the real 
value since, for a given sample, any source of unknown 
random error will tend to increase the measured velocity 
dispersion. The above value is significantly smaller than 
cTp = 2.2 ± 0.4 kms^^ found by Cote et al. 

Note that even a very crude unweighted computa- 
tion of the systemic velocity and projected velocity dis- 
persion will not produce significantly different results: 
namely, for the same eight stars, (14) = 25.95 kms~^ and 
CTp = 1.04 kms~^. In general terms, we can say that what- 
ever the sophistication of the weighting process used in the 
above computation, the resulting velocity dispersion will 
always be close to 1 kms^^. This is a robust result. 

What is the mass-to-light ratio corresponding to this 
new velocity dispersion dp value ? If only the cTp value 
changes, everything else being equal (King-Michie model, 
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Table 3. Membership probabilities from proper motions and positions in the CMD for the 11 stars with radial velocities 
implying membership of Palomar 13. 



ID QJ2000 5j2ooo At(a) At(<5) P(At) V, FLAMES U B V CMD 

[°"'] [masyr-i] [masyr"!] [%] [kms'i] [mag] [mag] [mag] 
(1) (2) (3) (4) (5) (6) (7) ^) (10) (11) 



38 


23 06 41.55 


12 46 58.5 


0.32 ±0.11 


0.05 ±0.19 


48 


112 


23 06 39.37 


12 47 35.1 


-0.31 ±0.16 


0.59 ±0.29 


89 


22 


23 06 37.31 


12 49 26.8 


-0.22 ± 0.07 


-0.02 ±0.10 


96 


24 


23 06 50.11 


12 47 14.9 


0.24 ± 0.19 


1.6 ±0.16 





41 


23 06 49.95 


12 45 28.7 


-0.05 ± 0.12 


-0.25 ± 0.12 


79 


79 


23 06 44.30 


12 45 49.4 


0.22 ±0.16 


-0.36 ±0.18 


65 


188 


23 06 45.76 


12 46 41.4 


-0.18 ±0.20 


0.09 ±0.23 


99 


37 


23 06 48.51 


12 46 19.8 


0.00 ± 0.09 


0.32 ± 0.09 


88 


127 


23 06 31.29 


12 43 36.0 


-4.4 ± 0.77 


-6.2 ± 0.54 





66 


23 06 48.28 


12 45 47.2 


-0.41 ±0.13 


0.03 ±0.14 


88 


88 


23 06 43.82 


12 45 36.4 


-0.30 ±0.18 


-0.20 ±0.20 


93 



core radius, and central surface brightness, see Table 3 
from Cote ct al.), the valnc of the M/Ly ratio varies as 
the squared ratio of the old to new velocity dispersions. 
Consequently, the new M. / Ly is equal to the old M. jLy 
divided by (2.2/0.92)^ = 5.8, namely, MjLv = 40/5.8 ~ 
7, a value significantly closer to the usual value M. jLy ~ 3 
for globular clusters. 

Since wc have a rather small sample, we explore briefly 
the sensitivity of the above results to various subsamples. 

4.2. Sample of five stars with two FLAMES K 
measurements 

Similar arithmetics applied to the five Palomar 13 mem- 
ber stars with at least two FLAMES measurements gives 
the values (y^) = 25.40 ± 0.32 kms'^ and cxps = 0.60 ± 
0.27 kms~^. This again may be considered as an upper 
limit to the real value. From (2.2/0.60)^ = 13.4, the new 
M./Ly = 40/13.4 ~ 3.0, a value identical to the usual 
value for globular clusters. 

4.3. Sample of six stars witli FLAMES and Kecl< K 
measurements 

Six of the above eight Palomar 13 stars have Vr measure- 
ments from both FLAMES and Keck studies. 

The six Vr measurements from this FLAMES study 
bring the values (Vrep) = 26.14 ±0.44 kms"'^ and (Tpep = 
0.99±0.34 kms~^. The corresponding A4/Ly would scale 
by (2.2/0.99)2 = 4.9, namely, M/Ly = 40/4.9 ~ 8. 

The six measurements from the Keck study bring 
the values {VrSK) = 24.60±1.16 kms'^ and c7p6K = 2.79± 
0.83 kms~^. The corresponding M/Ly would scale by 
(2.2/2.79)2 = 0.62, namely M/Ly = 40/0.62 ~ 60. 

For the same stars, FLAMES K values produce a 
smaller velocity dispersion than the Keck values, indi- 
cation of a possible better stability and better zero-point 
calibration for the simultaneous FLAMES measurements 
than for the independent Keck ones. 



21.39 18.432 ±0.020 18.335 ± 0.010 17.808 ± 0.007 RRL 

24.36 20.134 ±0.042 20.070 ± 0.021 19.276 ±0.014 SGB 

24.82 18.148 ±0.023 17.761 ±0.012 16.897 ±0.007 RGB 

25.39 18.334 ±0.020 17.983 ±0.011 17.005 ± 0.007 off-r? 

25.76 18.854 ±0.022 18.658 ±0.011 17.795 ±0.007 RGB 

25.76 19.418 ±0.027 19.386 ±0.014 19.021 ±0.012 BS? 

25.96 20.535 ±0.058 20.529 ± 0.024 19.834 ± 0.018 SGB 
26.32 18.757 ±0.023 18.537 ±0.011 17.625 ± 0.007 RGB 
26.44 21.740 ±0.185 20.579 ±0.034 19.305 ± 0.018 off-r 
26.69 19.427 ±0.028 19.361 ±0.014 18.583 ± 0.010 SGB 

27.97 19.842 ±0.033 19.727 ±0.016 18.970 ±0.011 SGB 



4.4. Overlap of Keck and FLAMES sample fields 

Because of the positioning of FLAMES fibres, which pre- 
vents the measurements of stars with small separations, 
the star sample acquired during this commissioning ob- 
serving run is not optimal for the study of Palomar 13. 
We do not have numerous member stars within the 
tidal radius of the cluster as determined by Cote et al. 
Following our three criteria (radial velocities, proper mo- 
tions, and CMD photometry), there seems to be members 
of Palomar 13 well outside the nominal tidal radius, as 
already noticed by Cote et al. 

Could our small velocity dispersion value originate 
from the fact that we have stars on average more distant 
from the cluster center than those in the Keck study ? This 
is not the case. We do not observe any central increase in 
velocity dispersion in Keck data since the stars contribut- 
ing strongly to the velocity dispersion lie all outside the 
radius 60 ". 

5. Final remarks 

Although small, the present sample of radial velocities 
measured with FLAMES provides velocity dispersion and 
mass-to-light ratio estimates which diff'er significantly 
from those obtained from Keck radial velocities. 

In both, FLAMES and Keck samples, the detection of 
binary stars is difficult. Cote et al. have good time cover- 
age (up to four epoch) but few stars measured more than 
once, while FLAMES has two measurements for most stars 
but no time coverage. Two stars, #66-41 and #37-72 are 
responsible for significant differences between FLAMES 
and Keck results. Though these objects have not shown 
significant variability in, respectively, three and four Keck 
measurements and fall in the Vr range of the cluster in 
FLAMES data, we must accept them as binary candi- 
dates. Based on the proportion of blue stragglers in CMD 
of Palomar 13, approximately 25% of evolved stars are bi- 
nary systems, a figure well compatible with two binaries 
found and simulations carried out by Cote et al. 
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The difference between FLAMES and Keck results 
casts doubt on the previous Keck velocity dispersion de- 
termination. The most simple and likely explanation of 
the very high Keck A4 / Ly ratio is the Keck sample con- 
tamination by binaries since the determination of velocity 
dispersion even from the larger (still small) sample from 
Keck is prone to be biased by the presence of few undected 
binaries. 

A complete and more detailed analysis of these 
FLAMES radial velocities along with the old and some 
new radial velocities from Keck is under preparation, with 
extensive simulations of such observations. A careful at- 
tention will be given to the surface density profile of 
Palomar 13 which anyway appears to be anomalous among 
Galactic globular clusters. Is there anything genuinely pe- 
culiar with Palomar 13 ? A significant effort (new observa- 
tions and new analysis) is needed to answer the following 
question: Is Palomar 13 the host of a large fraction of bi- 
naries which may infiate its velocity dispersion ? With a 
good observing strategy, FLAMES is currently the best 
instrument to answer this question. 
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his work in the development of the GIRAFFE Data Reduction 
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